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We report analytic solutions of a recently discovered quasi-exactly solvable model consisting of 
two electrons, interacting via a Coulomb potential, but restricted to remain on the surface of a 
D-dimensional sphere. Polynomial solutions are found for the ground state, and for some higher 
[L < 3) states. Kato cusp conditions and interdimensional degeneracies are discussed. 
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I. 



INTRODUCTION 



A quasi-exactly solvable model is one for which it is 
possible to solve the Schrodinger equation exactly for a 
finite portion of the energy spectrum [1]. In quantum 
chemistry, a famous example of this is the Hooke's law 
atom [2 5, which consists of a pair of electrons, repelling 
Coulombically but trapped in a harmonic external po- 
tential. This model and others [6HT2] have been used ex- 
tensively to test various approximations |13H20j within 
density functional theory (DFT) [2"TH2"3"] and explicitly 
correlated methods [2"lTf2"5] . 

We have recently discovered [29] that a pair electrons, 
repelling Coulombically but constrained to remain on 
the surface of a 2?-sphere of radius R jTTJ I5UH5T] is also 
quasi-exactly solvable and we have called this system T>- 
spherium. (We adopt the convention that a 2?-sphere is 
the surface of a (2?+ 1) -dimensional ball.) We have shown 
that the Schrodinger equation for the 1 5 l and the 3 P 
states of 2?-spherium can be solved exactly for a count- 
ably infinite set of R values and that the resulting wave 
functions are polynomials in the interelectronic distance 
< u = |n - r 2 \ < 2R. 

In this article, we extend our earlier results [29] to 
higher angular momentum (up to L = 3) states of T>- 
spherium (T> > 2) for both the singlet and triplet mani- 
folds. The T> = 1 case is anomalous and, for brevity, is 
not discussed here. We use atomic units throughout. 



where H, \ and \1/ are the spin, angular and interparticle 
wave functions, respectively, and Si and S~2^ are the spin 
and hyperspherical coordinates [33] of the i-th electron. 
The singlet and triplet wave functions are given by the 
familiar HO] forms 



1 -(si,s 2 



1 



[a(s 1 )/3(s 2 ) - /3(si)a(s 2 )} 



v/2 

{a(si)a(s 2 ), 
^[a(«i)^(sa) + /9(*i)a(s a )], 
/3(si)/3(s 2 ). 

The angular part is associated with an energy 

A 



(3) 
(4) 



A 



v - 1) 



V - 1) 



(5) 
(6) 



II. 



WAVE FUNCTION 



where t\ and £ 2 are the angular momentum quantum 
numbers of the corresponding one-electron configuration 
(s = 0, p = 1, d = 2, / = 3, . . . ). The functions Xi 
which are dependent on the nature of the state considered 
[38l [41] . are gathered in Table |TJ where 0i £ [0,tt] and 
4>i E [0, 27r] are the (T> — l)-th and 2?-th hyperspherical 
angles of the electron i. The corresponding one-electron 
configurations are also reported. In Table[l] the suffixes e 
(even) and o (odd) are related to the parity of the states, 
which is given by (— l) ll+l2 . Hence, we label the states 
with the notation 13 L C '°, where L = S,P,D,F, . . . 



The Hamiltonian of 2?-spherium is 

X = -lW + Vl) + l (1) 

where the two first terms represent the kinetic contribu- 
tion of each electron, and w -1 is the Coulomb operator. 

Following Breit [38], we write the total wave function 
as the product 

s 2 }, {n 1} n 2 }, «) = h(si, s 2 ) x (n 1) n 2 )*(«), (2) 
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III. POLYNOMIAL SOLUTIONS 

Substituting the ansatz ([2| into the Hamiltonian 
yields the Schrodinger equation 

/ u 2 U 2 * / Su 1\# * _ 

V4R2 ~ J + [aR* ~ ~^a) + u ~ ' ( ) 

where the parameters S and 7 are tabulated for each state 
in Table U 

By introducing the dimensionless variable x = u/2R, 
Eq. (0 can be recast as a Heun's differential equation 
[4"2"] with singular points at x = — 1,0, +1. Following 
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TABLE I. Ground state and excited states of D-spherium 



State 


Configuration 


X(fil,fi 2 ) 


r 




7' 


-1 


A 


K 


Transformation 


1 ftC 

A 


2 

S 


1 


2D — 1 


D - 


- 1 








3 pc 




sp 


COS 6\ — COS 02 


2D + 1 


D - 


r 1 


D/2 


1 




1 po 


sp 


COS 6i + COS 02 


2D + 1 


D - 


- 1 


D/2 





3 D° 


3 pc 


P 2 


sin 6i sin 62 sin ((j>i — $2) 


2D + 3 


D - 


h 1 


D 


1 




3 D e 


sd 


X3po • xipo 


2D + 3 


D- 


r 1 


D + 1 


1 


1 pe 


l D o 


pd 


X3po ■ X3 P o 


2D + 5 


D- 


h3 


3D/2 + 1 


2 




3 D o 


pd 


XI po ■ X3po 


2D + 5 


D- 


r 1 


3D/2 + 1 


1 




1 pe 


Pf 


X3pc • X 3 D° 


2D + 7 


D- 


h3 


2D + 3 


2 





the known solutions of this equation [43] , we seek wave 
functions of the form 



*(«) = £ 



a k u k , 



(8) 



fe=0 



and substitution into Q yields the three-step recurrence 
relation 



7 



(fc + 2) [(fc+ 1)7 + 1] 



AE 2 



E 



(9) 



with the starting values a = 1 and a! = 7. 

If the series Q terminates at a certain k 
obtain the exact wave function 

n 
k=Q 



n, we 



(10) 



for a particular radius i?„ !m and energy E n ^ m . This is 



an nth degree polynomial with m nodes between and 
2R (0 < m < [ ,J ^\) and requires that a„+i and a n+ 2 
vanish. If a n +\ = 0, the relation 

s-r 



(11) 



ensures that a„+2 = 0. For a given n, the energies 
are thus determined by finding the roots of the equation 
a n+ i = 0, which is a polynomial in E, of degree [ n Y^\ ■ 

For the 1 D C state, we have not been able to obtain 
polynomial solutions because the Hamiltonian ([!]) is not 
separable using the ansatz ^ and the wave function sat- 
isfies exchange coupled equations [44]. This applies also 
to some other states of higher angular momentum. 
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FIG. 



1. Energy of the S, P and D states of 3-spherium. 

< Spa < Ipo < 3pc < 3 D c < l D o < 3 D oy quag j_ 



exact solutions are shown by the markers. 



in Ref. [29] . Numerical values of the energies and radii 
for the higher angular momentum states can be deter- 
mined using the interdimensional degeneracies (see Sec. 



IVE) 



For any given state, as n increases, the radius increases 
and the energy decreases. The opposite behavior is ob- 
served with respect to m. Futhermore, as R (or, equiv- 
alently, n) increases, the electrons tend to localize on 
opposite sides of the sphere due to the dominance of the 
Coulomb interaction as the density decreases [Ml 155] , 
Such Wigner crystallization [45] has also been observed 
in other systems [3 [101 [46] . 

The energies of the S, P and D states (m = 0) for 
3-spherium are plotted in Fig. [T] (the quasi-exact solu- 
tions are indicated by markers), while density plots of 
2-spherium (n = 1 and m = 0) are represented on Fig. 
® 



IV. 



RESULTS AND DISCUSSION 



Numerical values of the energies and radii for t he 1 P° 
and 3 P C states are reported in Tables [Tl] and III Tables 
containing results for the 1 S' e and 3 P° states can be found 



A. Natural/unnatural parity 

In attempting to explain Hund's rules [27] and the "al- 
ternating" rule [1H1 HO] (see also [SO] [ST]), Morgan and 



3 



TABLE II. Radii R n ,m and energies £7 n>m for 1 P° states of two electrons on a D-sphere (2?=2,3,4) 







= 2 






2? = 3 






V = 4 




n/m 





1 2 


3 





1 2 


3 





1 2 


3 





1 


1.118 








1.871 








2.598 










2 


3.162 








4.637 








6.083 










3 


6.226 


1.656 






8.376 


2.520 






10.52 


3.303 






Radiu 


4 


10.30 


4.232 






13.11 


5.888 






15.93 


7.440 






5 


15.38 


7.847 


2.159 




18.84 


10.21 


3.127 




22.32 


12.49 


3.966 






6 


21.46 


12.49 


5.246 




25.57 


15.53 


7.077 




29.72 


18.51 


8.732 






7 


28.54 


18.14 


9.397 


2.639 


33.30 


21.84 


11.97 


3.707 


38.11 


25.51 


14.39 


4.599 




8 


36.63 


24.80 


14.59 


6.222 


42.03 


29.15 


17.86 


8.220 


47.49 


33.49 


21.01 


9.976 




1 


1.000 








0.5000 








0.3333 










2 


0.3000 








0.1860 








0.1351 








>> 


3 


0.1355 


1.914 






0.09622 


1.063 






0.07460 


0.7562 






bX 


4 


0.07541 


0.4467 






0.05820 


0.2884 






0.04731 


0.2168 






n 


5 


0.04757 


0.1827 


2.414 




0.03874 


0.1319 


1.406 




0.03260 


0.1041 


1.033 






6 


0.03257 


0.09620 


0.5450 




0.02753 


0.07468 


0.3594 




0.02378 


0.06129 


0.2754 






7 


0.02363 


0.05851 


0.2180 


2.764 


0.02051 


0.04771 


0.1587 


1.656 


0.01808 


0.04035 


0.1267 


1.241 




8 


0.01789 


0.03903 


0.1127 


0.6200 


0.01585 


0.03295 


0.08781 


0.4144 


0.01419 


0.02853 


0.07252 


0.3215 



TABLE III. Radii R„ >m and energies E„ iTn for 3 P C states of two electrons on a D-sphere (2?=2,3,4) 









V = 


2 






V = 


3 






V = 


4 






n/m 





1 


2 


3 





1 


2 


3 





1 


2 


3 




1 


2.291 








3.000 








3.708 










2 


5.477 








6.892 








8.307 








CO 


3 


9.616 


3.006 






11.72 


3.748 






13.84 


4.478 






3 

t3 


4 


14.73 


6.851 






17.52 


8.334 






20.32 


9.797 






p3 


5 


20.84 


11.62 


3.676 




24.30 


13.82 


4.453 




27.78 


16.01 


5.208 






6 


27.94 


17.35 


8.156 




32.07 


20.26 


9.708 




36.22 


23.15 


11.22 






7 


36.04 


24.05 


13.54 


4.315 


40.83 


27.66 


15.84 


5.128 


45.65 


31.25 


18.10 


5.909 




8 


45.13 


31.75 


19.86 


9.412 


50.58 


36.04 


22.90 


11.03 


56.07 


40.32 


25.89 


12.60 




1 


0.3333 








0.2500 








0.2000 










2 


0.1333 








0.1053 








0.08696 








>. 


3 


0.07300 


0.7472 






0.06002 


0.5874 






0.05093 


0.4862 









4 


0.04607 


0.2131 






0.03908 


0.1728 






0.03390 


0.1459 






a 


5 


0.03166 


0.1019 


1.018 




0.02752 


0.08503 


0.8196 




0.02429 


0.07315 


0.6913 






6 


0.02306 


0.05983 


0.2706 




0.02042 


0.05117 


0.2228 




0.01829 


0.04477 


0.1906 






7 


0.01752 


0.03932 


0.1242 


1.222 


0.01575 


0.03432 


0.1046 


0.9983 


0.01428 


0.03046 


0.09079 


0.8522 




8 


0.01375 


0.02778 


0.07096 


0.3161 


0.01251 


0.02464 


0.06102 


0.2629 


0.01145 


0.02215 


0.05370 


0.2269 



Kutzelnigg [52"H5"i] have proposed that the two-electron 
atomic states be classified thus: a two-electron state, 
composed of one-electron spatial orbitals with individual 
parities (— l)^ 1 and (— l)^ 1 and hence with overall pari- 
ties (— lY 1+l2 , is said to have natural parity if its parity 
is (— 1) L . [■ ■ ■] If the parity of the two-electron state is 
— (— 1) L , the state is said to be of unnatural parity. |54j 

After introducing spin, three classes emerge. In a 3- 
dimensional space, the states with a cusp value of 1/2 [53J 



ISrJ] are known as the natural parity singlet states, those 
with a cusp value of 1/4 [57] are the natural and unnatural 
parity triplet states, and those with a cusp value of 1/6 
[5 2) , are the unnatural parity singlet states. 

In previous work [S5], we have observed that the 1 S C 
ground state and the first excited 3 P° state of 3-spherium 
possess the same singlet (1/2) and triplet (1/4) cusp con- 
ditions as those for electrons moving in three-dimensional 
physical space and we have therefore argued that 3- 
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FIG. 2. Density plots of the S, P and D states of 2-spherium. 
The squares of the wave functions when one electron is fixed at 
the north pole are represented. The radii are %/3/2, VT5/2, 
y/E/2, V21/2, V21/2, 3/5/2 and 3^3/2 for the 3 P°, 
1 P°, 3 P C , 3 L» C , 1 D° and 3 D° states, respectively. 



spherium may be the most appropriate model for study- 
ing "real" atomic or molecular systems. This is sup- 
ported by the similarity of the correlation energy E c of 
3-spherium to that in other two-electron systems. In- 
deed, it can be shown that, as R — > 0, E c approaches 
—0.0476, which is close to the corresponding values for 
the helium-like ions (—0.0467) [58], the Hooke's law atom 
(-0.0497) [5S], and two electrons in a ball (-0.0552) 06J. 

Most of the higher angular momentum states of 3- 
spherium, possess the "normal" cusp values of 1/2 and 
1/4. However, the unnatural 1 D° and 1 F C states have 
the cusp value of 1/6. 



B. First-order cusp condition 

The wave function, radius and energy of the lowest 
states are given by 

¥i,o(i0 = l + 7u, R la = 7-> E ifi=l, (12) 

47 

which are closely related to the Kato cusp condition [56] 
*'(0) 



tt(0) 



7- 



(13) 



We now generalize the Morgan-Kutzelnigg classifica- 
tion [53] to a 2?-dimensional space. Writing the interpar- 



ticle wave function as 
*(tt) = H 

we have 



2k + V - 1 



0(u 2 



(14) 



K = for natural parity singlet states, 

K = 1 for triplet states, (15) 

k — 2 for unnatural parity singlet states. 

The labels of the P-spherium states are given in Table [T] 
C. Second-order cusp condition 



The second solution is associated with 

(6 + 2) 



*2,o(«) = *l,o(«! 



R 2,0 = 



^2,0 — 



2 7 (5 + 2) +45- 
( 7 + 2)(<5 + 2)-l 



2 7 
7(5+1) 



( 7 + 2)(* + 2)-l 
For I?-spherium, the second-order cusp condition is 
*"(0) 1 ( 1 



*(0) 



22? V V - 1 



E 



(16) 
(17) 
(18) 

(19) 



Following ( 19 ), the classification ( 15 ) can be extended to 



the second-order coalescence condition, where the wave 
function (correct up to second-order in u) is 



*(u) 



1 



2k 
,,2 



V-l 



1 



2(2k + V) \2k + V- 1 
Thus, we have, for T> = 3, 

*"((). = I \ { f~ E ^ 

1 10 u 



— E \ + 0(u ). (20) 



*(0) 



14 V 6 



for k = 0, 
for k = 1, 
for k = 2. 



(21) 



For the natural parity singlet states (k = 0), the 
second-order cusp condition of 3-spherium is precisely 
the second-order coalescence condition derived by Tew 
[60) , reiterating that 3-spherium is an appropriate model 
for normal physical systems. 



D. Third-order cusp condition 

The third-order cusp condition of 3-spherium is 
<F'"(0) 1 ( \ 



*(0) 



18 V 8 



£+ l^ 



(22) 



which is similar, but not strictly equivalent, to the one 
derived by Tew [BD], due to the Independence of (22 1. 



The generalization to T> dimensions is straightforward. 
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E. Interdimensional degeneracies 



As shown in Table [TJ many states of 2?-spherium have 
the same Hamiltonian |7]) as lower angular momentum 
states of (V + 2)-spherium. 



Using the transformation (T>, L) — > (D + 2, L — 1) (see 
Table [M), one can see that the Hamiltonian of the 3 P C , 
1 D°, 3 D° and 1 i ?c states for a given value of I? are respec- 
tively identical to those for 1 5°, 3 P°, 1 P°, and 3 D° states 
at O + 2. The transformation (D, L) (V + 2, L - 1), 
preserves the parity of the states, but "flips" the spin 
configuration, thereby increasing by one unit the value 
of k. In 2?-spherium, we note that the Hamiltonians of 
the 3 P C and 3 D C states are identical. 



Similar interdimensional degeneracies, first noticed by 
van Vleck [5T], have been observed for various systems 

BHE3HB5]. 



V. CONCLUSION 

In this article, we have reported exact solutions of a 
Coulomb correlation problem, consisting of two electrons 
on a 2?-dimensional sphere. The Coulomb problem can 
be solved exactly for an infinite set of values of the radius 
R for both the ground and excited states, on both the 
singlet and triplet manifolds. The corresponding exact 
solutions are polynomials in the interelectronic distance 
u. 

The cusp conditions (up to third-order in the inter- 
electronic distance), which are related to the behavior 
of the wave function at the electron-electron coalescence 
point, have been analyzed and classified according to the 
natural or unnatural parity of the state considered. 

Finally, we have shown seen that, as in other one-, two- 
or three-electron systems, there exist interdimensional 
degeneracies between some of the states of P-spherium. 
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